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Abstract 

We discuss the experimental requirements for a mass hierarchy measurement for 13 = 
using muon neutrino disappearance. We find that a specially optimized neutrino factory at 
L ~ 6 000 km could do this measurement using extreme luminosities. In particular, we do 
not require charge identification for this purpose. In order to measure the mass hierarchy for 
more adequate luminosities, we explore the capabilities of low energy narrow band off-axis 
beams, which have relatively more events at low energies. We find that, in this case, the 
energy resolution of the detector quickly becomes the limiting factor of the measurement, 
and significantly affects the baseline optimization for determining the mass hierarchy. 
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1 Introduction 



Among the goals of next-generation neutrino experiments are a measurement of the magni- 
tude of the #13 angle of the recently revealed lepton mixing matrix (using the PDG param- 
eterization of the mixing matrix [1]) and the determination of the so-called neutrino mass 
hierarchy. 

Determining the neutrino mass hierarchy is equivalent to establishing how the neutrino mass 
eigenvalues, properly defined, are ordered. It is customary to define the neutrino masses- 
squared mf, m|,m| such that m\ < m\ and Am^ < lAmJ-J, where Am| = mf — m|. The 
case tttJ > m\ (Am^ > 0) is referred to as a normal mass hierarchy, while the other logical 
possibility, m\ < m\ (Am^ < 0), is referred to as an inverted mass hierarchy. Current data 
on neutrino oscillations have allowed the determination of Am^ and jAm^l, but not the 
sign of Am^. 

In the case of "large" #13 values, several neutrino oscillation studies have been performed 
in order to address the optimal means for addressing the neutrino mass hierarchy. These 
include studies of —>■ v e transitions in matter (see e.g. Refs [2-12]), v e — > transitions 
in matter (see e.g. Refs. [13-21]), precision studies of the neutrino flux from supernovae 
(see e.g. Refs. [22-26]), disappearance of atmospheric neutrinos [27-29], and combined 
analyses of and v e disappearance in vacuum [30,31]. In the case of "small" 6*13 values, on 
the other hand, the issue of determining the neutrino mass hierarchy has been neglected, 
until recently. 

In [30], it was shown that in the limit #13 — > 0, the neutrino mass hierarchy can be deter- 
mined via precision studies of disappearance under rather extraordinary conditions. The 
discussion in [30] was, however, purely theoretical, and no attempt was made to determine 
whether such studies could be performed in practice. Here, we discuss the experimental re- 
quirements for establishing the neutrino mass ordering in the limit #13 — > via long-baseline 
studies of (and u^) disappearance. 

First, we revisit in Sec. |5]the analysis performed in [30] in order to discuss the conditions that 
need to be met in order to determine the neutrino mass hierarchy via neutrino oscillations. 
We review that it is imperative to combine precise measurements of the — > survival 
rate at different values of L/E and that at least one of the L/E values should satisfy 
A17121L/E ~ 1. In Sec. El we study in more detail what turns out to be the issue that more 
significantly limits our ability to determine the neutrino mass hierarchy — the detector 
energy resolution in the case of low-energy (hundreds of MeV) muon-type neutrinos. 

In Sec. we study the ability of "conventional" neutrino beams and detectors to study the 
neutrino mass hierarchy via disappearance. After convincing ourselves that it is very 
unlikely that the neutrino mass hierarchy will be uncovered with these types of beams and 
detectors, we discuss, in Sec. El how the situation could change with improved detector 
energy resolution. In Sec. we summarize our results and conclude. 
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2 Physical mechanism 



In the limit #13 — ► 0, which we assume applies henceforth, oscillations involving electron-type 
neutrinos (or antineutrinos) in the production or detection stage are sensitive only to one 
neutrino mass-squared difference: Ara^. Therefore, they cannot convey any information 
about the sign of Am^, and, in order to learn about the neutrino mass hierarchy, one is 
forced to turn to the oscillation of muon-type neutrinos and antineutrinos. 

In a nutshell, P^ 1 in vacuum is sensitive to the mass hierarchy in the following way (for 
many more details, see [30,31]). We can express 

where A,B,C are functions of the lepton mixing angles and Ajj = Am|-/2S. In the case 
of a normal hierarchy, |A 31 | > |A 32 |, while in the case of an inverted one the converse is 
true. Hence, if (a) one can determine that three different frequencies contribute to P^ and 
(b) B 7^ C (and both are known), the mass hierarchy is trivially determined. In the case 
of a normal mass hierarchy, the amplitude associated with the largest |Ay| is B, while in 
the case of an inverted mass hierarchy the amplitude associated with the largest |Ajj| is 
C. Condition (a) summarizes the first challenge one needs to meet: the experimental setup 
should be such that Am^ and Am§ 2 effects are distinguishable. 

In practice, the situation is more involved. The reason is that [Am^l is measured at the 
same time that the neutrino mass hierarchy is determined. It turns out that [30] , for a fixed 
value of L/E, for every Am^ = Am^ > 0, there is at least one different Am^ = Am^ < 
such that 

P w (Am^ = Am 2 +) = P^(Am 2 31 = Am^). (2) 
Furthermore, if A 21 L <C 1, 

Am 2 3 i = -Am^ + 2Am 2 21 cos 2 6 12 , (3) 

independent of L and E. The results above do not depend on whether the neutrinos are 
oscillating in vacuum or constant matter densities. 

Throughout, we use the following neutrino oscillation parameters (c/., Refs. [32-36]), unless 
stated otherwise: 

sin 2 (2#i 2 ) = 0.83, Am 2 21 = 8.2 10~ 5 eV 2 , 

sin 2 (2# 13 ) = 0, sin 2 (2# 23 ) = 1.0, (4) 

Note that, since we assume that 6*i 3 is exactly zero, 5qp is not a physical observable. Fig. [T] 
depicts P w as a function of E for L = 295 km, the oscillation parameter values in Eq. (J3J), 
and three different values of Am 2 3 : 2 Am 2 x = +2.2 10~ 3 eV 2 (solid dark [black] curve), 

■""We denote P a @ as the probability that a neutrino with energy E produced in a weak eigenstate v a is 
detected, having propagated a distance L, as a weak eigenstate up. P & p is the equivalent for antineutrinos. 

2 According to recent analyses of atmospheric data [35], lAmfJ is constrained to be, at the 90% level, 
between 1.9 10" 3 eV 2 and 3.4 1(T 3 eV 2 , while the best fit point is located between 2.0-2.5 10" 3 eV 2 . 
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Figure 1: P w as a function of E for L = 295 km, the oscillation parameter values in 
Eq. and three different values of Am 2 3 : Am| x = 2.2 10~ 3 eV 2 (solid dark [black] curve), 
Am^ = -2.2 10 _3 eV 2 (solid light [green] curve), and Am^ = -2.08 10~ 3 eV 2 (dashed light 
[green] curve). 

Amj! = -2.2 l(T 3 eV 2 (solid light [green] curve), and Am^ = -2.08 10" 3 eV 2 (dashed light 
[green] curve). One can clearly see that there are two virtually identical curves: One asso- 
ciated to a normal mass hierarchy, the other to an inverted one. Furthermore, as discussed 
above, the two distinct jAm^l values that yield the same oscillation probabilities are shifted 
by the amount 2Am 2 , 1 cos 2 9\2 = 1.1610~ 4 eV 2 . A measurement of Am^j performed under 
these circumstances would yield, regardless of its precision, two degenerate values — one 
associated to a normal mass hierarchy, the other associated with an inverted one. It is in- 
structive to note that, if, for some magical reason, jAm^l were independently known with 
great precision, it would be "straight forward" to tell an inverted from a normal hierarchy 
by measuring P MM at L = 295 km (as can see by comparing the solid lines in Fig. 

In order to break the degeneracy, one must also study P^ at a different value of L/E. 
Furthermore, it is imperative that the value be "large," such that A21L is of order one, and 
Eq. does not apply. Numerically, 

/ Am 2 \ /0.6GeV\ f L \ 

Hence, A 2 iL > 1 implies long distances [L > 3 000 km) and small neutrino energies, 
(P<lGeV). 

Fig. |21 depicts as a function of E for L = 6000 km, the oscillation parameters in Eq. (@J), 
and the two distinct values of Am^ that yielded the same P w functions at L = 295 km: 
Am^ = 2.210~ 3 eV 2 (solid [black] curve) and Am 2 31 = -2.08 10~ 3 eV 2 (dashed [green] 
curve). At large energies (small A 2 ±L), the situation is identical to the L = 295 km case 
(Fig. HJ) - - the ? w functions are virtually identical when one compares the best fits to 
the L = 295 km data provided by the two hypothesis concerning the mass hierarchy. The 
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Figure 2: P w as a function of E for L = 6000 km, the oscillation parameters in Eq. (JJJ), 
and Amjj = 2.210" 3 eV 2 (solid [black] curve) and Am^ = -2.0810~ 3 eV 2 (dashed [green] 
curve). 

situation is quite different at low energies (order one A 2 iL). There, the expected values of 
are quite distinct, especially at very low energies. Hence, while the first few oscillation 
maxima are positioned at the same energy value (e.g., in the figure these are occurring at 
around 10.5 GeV, 3.5 GeV, 2.1 GeV, . . . ) for both hierarchy hypothesis, higher oscillation 
maxima disagree. If these higher oscillation maxima were experimentally resolved, one 
should be able to uniquely determine the neutrino mass hierarchy. This is the issue on 
which we concentrate in the next section. It will turn out that nontrivial information about 
the neutrino mass hierarchy will arrive, for L = 6000 km, from neutrino energies E < 2 GeV 
(or so). Note that while the difference between the position of the minima of the two curves 
is very small for E = 1 or 2 GeV energies (around 10 MeV), the cumulative effect of 
measuring the position of several minima will prove to be important (cf. Sec. EJ). 

In principle, more information can be obtained if one also explores Pp,p, the muon-type 
antineutrino oscillation probability. In Fig. El we depict P^p as a function of E for L = 
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Figure 3: Ppp as a function of E for L = 6000 km, the oscillation parameters from Eq. (J3J, 
and Am^ = 2.210~ 3 eV 2 (solid [black] curve) and Am^ = -2.0810~ 3 eV 2 (dashed [green] 
curve). 

6000 km, the oscillation parameters from Eq. (J3J), and Am^ = 2.2 10"3 eV 2 (solid [black] 
curve and Am 2 31 = -2.08 10~ 3 eV 2 (dashed [ green] curve). As in Fig. El (top), the two 
candidate hypothesis regarding the neutrino mass hierarchy yield (in this case only very 
slightly) different survival probabilities as a function of energy. Most relevant here is the 
fact that the behavior of neutrinos and antineutrinos is distinct, which means that the "fake 
solution" patterns are shifted into different directions. Hence, if one were to combine studies 
of both polarities, one should be able to better discriminate the two mass hierarchies. The 
origin of this distinction is due to nontrivial Earth matter effects, which are present for 
these very long baselines. We emphasize that these are "solar" matter effects, that render 
the effective #™ atter and the effective oscillation frequencies A^ atter and A^ atter distinct for 
neutrinos and antineutrinos (for detailed expressions, see, for example, Refs. [30,37]). Note 
that this effect of Earth matter is conceptually different from the one used to probe the 
mass hierarchy via «-> z/ e -oscillations in the case of "large" #13 values. 

3 Experimental limitations 

It is, mainly, a detector's energy resolution that determines its ability to "see" oscillation 
maxima and minima. As described above, it is this ability that will determine whether 
an experiment can determine the neutrino mass hierarchy via muon-type neutrino (and/or 
antineutrino) disappearance. 

The technique used to determine the neutrino energy depends on the detector and the 
detection process. It is always the case, however, that the neutrino is observed via charged 
current scattering off some target X: iv^) + X —>■ fi T + X', and the incoming neutrino 
energy is reconstructed from the measured muon energy and, sometimes, the measured or 
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Figure 4: P w averaged over AE = 80 MeV energy bins, as a function of E for L = 
6000 km, the oscillation parameters from Eq. (0J), and Am| x = 2.210 _3 eV 2 (open [black] 
circles) and Am^ = — 2.0810~ 3 eV 2 (solid [green] asterisks). 

inferred X' energy. Water Cherenkov detectors, for example, detect neutrinos via charged- 
current neutrino scattering on nucleons, and the incoming neutrino energy is reconstructed 
by assuming that the scattering process is quasi-elastic. In calorimeter-like detectors (such 
as MINOS), on the other hand, the recoil hadronic energy is often also measured, and 
included in the reconstruction of the incoming neutrino energy. 

As we are interested in the best achievable energy resolution for low-energy neutrinos, we 
emphasize here that, when experiments such as NOvA talk about effective energy resolutions 
AE/GeV ~ 3% — 6% x ^/P/GeV, they refer to the relevant analysis range where their main 
signal is (in the case of NOvA, around a couple of GeV). This estimate would translate, for 
a 1 GeV neutrino, into a 30 — 60 MeV energy resolution, and is dominated by uncertainties 
related to measuring the muon and recoil hadronic energies, understanding the physics of 
neutrino nucleon scattering in the nuclear environment, etc. However, there is a "physics" 
lower limit on the energy resolution as E decreases, dominated by the Fermi motion of 
nucleons inside the nuclei. For example, the energy resolution of the QE events in a water 
Cherenkov detector is already dominated by this Fermi motion, which leads to an effective 
energy resolution of about 85 MeV [38]. It is in this lower limit we are ultimately interested 
here, and it defines what we mean by "energy resolution": It is the minimum achievable 
energy resolution for sub-GeV neutrinos. 

Fig. |U depicts P MAt averaged over constant AE = 80 MeV energy bins (corresponding 
to an energy resolution of about ±40 MeV) as a function of E for L = 6000 km, the 
oscillation parameters from Eq. and Am^ = 2.2 lO- 3 eV 2 (solid [black] curve) and 
Am 2 x = -2.08 10~ 3 eV 2 (dashed [ green] curve). The averaging over energy bins is meant to 
illustrate the effect of a finite energy resolution: as can be seen in the figure, all "wiggles" 
below 700 MeV are erased, together with the distinction between the two different hierarchy 
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hypothesis (compare Figs. |Uand[2J). 3 

It is easy to estimate when the oscillations are "too fast" to be resolved. The difference 
between the energy of the nth and the n + 1th maxima (or minima) is 

E n — E n+ i 2 



E n 2n + 3' 

where, in vacuum, the nth maximum is defined by 4 



(6) 



lAmli \L , . 7r 



4 En < 2 " + 1 >r m 

For a fixed energy resolution AE, we can (from the requirement E n — E n+ i > AE) at best 
hope to probe n values smaller than 



Am\ x \L 
AnAE 



n max ~ \l 1 31 ' =8.2 



[Am|j \ / 80 MeV \ / L 
2.2 x 10- 3 eV 2 / \ AE ) V6000 km 



where we have assumed that n max 3> 1. It is clear that larger n values are accessible for 
larger values of L and smaller values of AE. Note that the maximum n max is only observable 
if the number of events at E nmgx (defined in Eq. (jZJ) is large enough and is located within the 
analysis range. This means that the nth maximum may in fact appear below the analysis 
range, in which case the sensitivity of the measurement is not predominantly determined 
by AE. Using Eq. © and Eq. ©, we have 



E^ocy/AmlAEL. (9) 

Thus, smaller values of AE correspond to smaller energies E nm ^. This is easy to understand 
since, for fixed L, the "faster" oscillations at lower energies can be resolved with a better 
energy resolution (cf., Fig. EJ. If, on the other hand, the baseline increases for fixed AE, 
n max will increase according to Eq. (jBJ) and will move to higher energies according to Eq. (JHJ) , 
because the oscillations become faster at the energy where it was before. In practice, 
statistics plays a crucial role in realistic experiments. For instance, if we assume a very 
narrow band beam which only has events close to E, then the constraint E ~ £„ max implies 
with Eq. (JSJ) that the optimum baseline L will scale as L oc (Am| x AE) -1 to resolve n max . 
On the other hand, the 1/L 2 drop of the flux and the shape of the spectrum might actually 
imply that it is more efficient to resolve n max — 1 at a shorter baseline (cf., Eq. (jHJ)). In the 
next two sections, we will see that it is a non-trivial interaction between these factors which 
determines the baseline optimization for narrow band beams. 

In summary, while the effect we are interested in is concentrated at the lowest energies, the 
finite detector energy resolution washes out the effect, and requires one to obtain information 
regarding the mass hierarchy from intermediate energies. In order to understand the con- 
sequences of the interplay between these two conflicting requirements, detailed simulations 
are required. These are described in detailed in the next sections. 

3 "Left-over" differences between the two hierarchy hypothesis at averaged out energies are an artifact of 
the fixed bin size of Fig. They are not observable in more realistic simulations, including those reported 
in the other sections of this paper. 

4 The distinction between oscillations in matter and in vacuum is, as far as these estimates are concerned, 
very small, and can be ignored. 



7 




Figure 5: The binned energy spectrum for the disappearances channels of our NB- 
Superbeam (left, L = 6 600km) and Neutrino Factory (right, L = 6 000km), where the 
energy resolution of the NB-Superbeam has been exaggerated (improved) by a factor of two 
to identify the wiggles. For the oscillation parameters, we use the values in Eq. (j3J), and 
Amfi = +2.2 10" 3 eV 2 . 

4 Conventional methods 



In this section, we focus on established beam and detector technologies, such as detector 
technologies currently used or in active preparation. In addition, we do not assume an energy 
resolution pushed to the extremes. We first discuss how one would design such experiments 
with respect to the mass hierarchy measurement, and then evaluate their performance. 

4.1 Description of experiments 

Following the results of the previous sections, we will concentrate on setups capable of 
measuring Amj^ via and/or z/ M disappearance with high precision (high luminosities) at 
(at least) two different values of L/E and at small energies E < 1 GeV. In particular, we 
are interested in the ability of such experiments to resolve high oscillation maxima, which 
means that they required very good energy resolution. We consider two possibilities for this 
measurement at two distinct values of L/E: 

1. Two Narrow Band (NB) beams at different values of L. 

2. One Broad Band (BB) beam with a wide range in E. 

In order to probe option (1), we assume a well-established detector and beam, namely, we 
combine an off-axis superbeam with a water Cherenkov detector, using the T2HK proposal 
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as an archetype. This experiment peaks at an energy of about 0.76 GeV, which is perfectly 
suited for our purposes. To study option (2), one could use either a broad band superbeam 
or a neutrino factory beam with much higher statistics. Here, we concentrate on a neutrino 
factory beam in order to test higher statistics signals. Note that, in principle, a broad band 
superbeam could replace two narrow band beams, but the detector requirements will be 
similar and the events in the relevant (low) energy range will rather decrease than increase. 
One potential candidate for such a broad band superbeam is the BNL-Homestake beam [39]. 

Narrow band superbeam and water Cherenkov detector 

For the narrow band superbeam (NB-Superbeam), we use a setup similar to the T2HK 
off-axis superbeam [38] with the experiment simulation from Ref. [17] using the GLoBES 
software [40]. This superbeam uses a target power of 4MW, a baseline of 295 km in its 
T2HK configuration, and a water Cherenkov detector with a fiducial mass of 1 Mton. For 
the running time and polarity, we choose eight years of neutrino running only unless stated 
otherwise. The energy resolution for the QE events is mainly determined by the Fermi- 
motion of the nucleons in the oxygen core, which means that one has an effective energy 
resolution of about <je — 85MeV. 5 Note that this energy resolution is approximately 
constant in energy. Compared to Ref. [17], we change the binning and sampling points 6 
in order to be able to resolve fast oscillations and higher oscillation maxima. We use 50 
(equally spaced) bins between 0.4 and 1.2 GeV and 75 (equally spaced) sampling points 
between 0.2 and 1.4 GeV to avoid aliasing effects from the energy resolution function. 7 This 
corresponds to a constant bin/sampling width of AE = 16MeV, which means that, in 
principle, energy resolutions up to AE = 8 MeV can be tested. However, this sampling 
width is not enough to calculate very fast oscillations at low energies. Therefore, we use the 
low-pass filter of GLoBES (c/., GLoBES manual) to average out fast oscillations already at 
the probability level, and average over the bin width AE = 16 MeV. Since we have now 
artificially increased our effective energy resolution by averaging at two different places, we 
subtract this probability-level averaging effect from the energy resolution oe such that we 
obtain an effective energy resolution a c g ~ ^/ (J e+ (16 MeV) 2 ~ 85 MeV (energy resolution 
type 2 in GLoBES). In Fig. |5] (left), we show a typical energy spectrum with an energy 
resolution a e g exaggerated by a factor of two (the bin-widths are 16 MeV, considerably 
smaller than the energy resolution). One can see that there are enough events/bin to contain 
significant statistical information, while the "wiggles" can be resolved for E > 700 MeV. 

5 For the experiment simulation, we use an energy resolution defined as a Gaussian map- 
ping from the incident neutrino energy E to the reconstructed neutrino energy E': R(E,E') = 

(o-(E)y/2ir)~ 1 exp ^— ^a^fJ) ) I^L where a(E) is the energy resolution AE as function of E. 

6 Number of points, where the oscillation probabilities are evaluated. 

7 Since the event rates at the lower and upper end of the spectrum are not exactly zero, events beyond 
the analysis range could contribute via the energy resolution function. Therefore, it is advisable to use a 
larger sampling range than bin range. 
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Wide band neutrino factory beam and TASD 

We choose a neutrino factory similar to Ref. [17] with 1.06- 10 21 useful muon decays per year 
and a total running time of four years in each polarity (corresponding to 5.3 ■ 10 20 useful 
muon decays per year and polarity for a simultaneous operation with both polarities). Since 
we want to have as many events as possible at low energies, we use a lower muon energy 
of = 15 GeV. Though the total event rate for this choice is much lower than for more 
"canonical" ~ 20 — 50 GeV, it turns out that the relative increase in the low energy events 
is favorable for the performance. For the detector, we use the Totally Active Scintillator 
Detector (TASD) from Ref. [21] with a fiducial mass of 50 kt, which is similar to the NOvA 
detector [41]. In principle, one could also use an iron calorimeter, but the energy resolution 
and efficiency of the TASD is much higher. Finally, we do not require charge identification, 
as it would drastically reduce the efficiencies at low energies. Instead, we effectively measure 
the disappearance rates = A MM + N Efi or A?? = + N eil , where the event rates N a/3 
correspond to P a $ convoluted with the corresponding beam spectrum, cross sections, etc. 
This implies that the detector cannot distinguish muons from antimuons. 

For the binning, we use a somewhat more sophisticated approach compared to the super- 
beam one described above, because the neutrino factory spans a very large energy range 
with oscillation probabilities changing slowly in the high energy range and very quickly in 
the low energy range. We divide the energy analysis range from 0.5 to 15 GeV into variable 
bins with sizes of 100 MeV until 3.5 GeV, 250 MeV until 9.5 GeV, and 500 MeV for the rest, 
where we include an additional bin over 15 GeV to reduce aliasing effects. For the muon en- 
ergy resolution, we use 0.03 ^/E/GeV. Since this corresponds to an energy resolution of only 
21 MeV at the lower end of the analysis range and the Fermi motion contribution is much 
higher, we add an estimated resolution of AE = 85 MeV via the filter feature of GLoBES 
(carbon should have a similar order Fermi motion contribution compared to oxygen). This 
yields an effective best energy resolution a cS ~ a/ (85 MeV) 2 + (16 MeV) 2 ~ 87 MeV at the 
lower end of the spectrum, and <T e ff ~ 100 MeV at the energy range where first significant 
contributions to the event rates come in. We show in Fig. |3] (right), a typical energy spec- 
trum for the neutrino factory. At low energies, the bins are chosen fine enough to resolve 
fast oscillation not averaged out by the energy resolution, but not too fine such that there 
are enough events to significantly resolve the left-most oscillation maxima. As for most 
broad band beams, the lower energies are suppressed by the beam spectrum. This implies 
that though the mass hierarchy effects become smaller for larger energies, the statistics 
rapidly increases and partially compensates for this decrease. Thus, compared to Fig. 
(left), the broad band beam's performance should be significantly determined by more than 
one oscillation maximum. 

4.2 Performance analysis 

Before we come to the analysis of the described experiments, we need to define our perfor- 
mance indicator: We define that we have sensitivity to the normal (inverted) mass hierarchy, 
if there is no solution with the inverted (normal) mass hierarchy fitting the original solution 
below the chosen confidence level. Therefore, the mass hierarchy sensitivity is in practice 
determined by the sgn(Am| 1 )-degeneracy [42]. Since it turns out that for the disappearance 
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Figure 6: Ax 2 for the test of the normal hierarchy as function of two baselines for the 
different experiment combinations as given in the plot labels. The different contour curves 
correspond to Ax 2 = 0.05 (dashed), 0.1 (solid), 0.15 (dashed), 0.2 (solid) etc., whereas the 
different shadings correspond to Ax 2 increased in steps of 0.1. The diagonal lines represent a 
single baseline only, i. e., all detector mass put to one baseline. For the oscillation parameters, 
the values in Eq. (JU) are used, and Am^ = +2.2 10~ 3 eV 2 . 

channel in the limit sin 2 (2# 13 ) = the position of this degeneracy is mainly determined by 
a different value of Am^ [30] , we locate it in Am^ -direction first and then marginalize over 
all oscillation parameters (including Am^) using this starting point. In the marginalization, 
we assume an external precision of 5% for Amfj and 10% for 9i%, which is a conservative 
choice for the time when the experiments are being analyzed (c/., e.g., Ref. [33]). Further- 
more, we include a matter density uncertainty of 5 % [43-45]. We do not include external 
constraints for sin 2 (2^13). 

Since we know that, in principle, we need at least two different values of L/E to resolve 
the neutrino mass hierarchy, we combine a priori two different experiments/baselines to 
demonstrate the effects. In Fig. we show the Ax 2 for the test of the normal hierarchy as 
function of two baselines for the different experiment combinations as given in the plot labels. 
We first note that the absolute Ax 2 in these figures is at maximum one, which means that 
the experiments considered are not sufficient for the measurements and require substantial 
luminosity upgrades. However, these figures are useful to discuss the optimization of these 
experiments for the mass hierarchy measurement. 

For the combination of two NB-Superbeam baselines (left panel of Fig. |)J), we find the 
optimal performance for the combination 295 km (lower end of plot) + 2 700 km. This is 
what one would expect from a narrow band beam: One baseline is not sufficient to resolve 
the mass hierarchy. Instead, one needs the precise Am^ measurement at two very distinctive 
values of L/E. It is not very surprising that the baseline L = 295 km is good for a precision 
measurement of Am^ because it corresponds to the original T2HK setup. This baseline 
basically determines the value of the "fake" Am| x of the inverted sign solution, which is 
then compared to the "fake" Am^ at the longer baseline where solar effects contribute. 
The inconsistency between these two Am^'s is then represented by the Ax 2 - As we will see 
later, the length of the long baseline depends on the energy resolution. It is determined by 
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Figure 7: Potential detector locations L = 6 000 km for three different potential neutrino 
factory laboratories. 

a balance between several competing factors. As discussed in Sec. the energy resolution 
(and the beam energy profile) determines the smallest energy where one can still "measure" 
Amf 3 . With that in mind, the optimum baseline is such that the largest possible oscillation 
maximum can be resolved with enough luminosity (larger n max prefers larger values of L, 
where more statistical significance is obtained than for shorter baselines; cf., Eq. (|5jl). 

As already mentioned above, the luminosity is still not at all sufficient for a significant 
result: In this case, a factor of 30 luminosity upgrade would be necessary for a 3a signal. 
This corresponds to a total luminosity of ~ 480Mtonyr (4MW beam power), which is not 
achievable in practice. A somewhat better result can be achieved by using neutrinos and 
antineutrinos (2 yr + 6 yr) at one baseline L ~ 2 700 km instead of two different baselines. 
In this case, a factor of 20 luminosity upgrade would be sufficient. In summary, two narrow 
band superbeams do not seem to have sufficient statistics for this measurement, but the 
theoretical effect is, in principle, observable. 

The ideal instrument to boost the statistics of such a measurement is a neutrino factory. 
Fig. |H1 (middle) depicts the reach of combining data from the NB-Superbeam and the neu- 
trino factory we introduced above. A neutrino factory at L ~ 6 000 km would clearly 
enhance the potential and complement the L ~ 2 700 km baseline of the NB-Superbeam. 
Note that the A% 2 of the combination of these two experiments is higher that the sum of 
the Ax 2 's of the neutrino factory alone and the superbeam combination alone, which means 
that there is real synergy between these two experiments beyond just adding their statistics. 
For two neutrino factories (case depicted in Fig. (right)), it is eventually best to put all 
detector mass at L ~ 6 000 km. As one may expect, the broad band nature of the beam 
itself has enough oscillation maxima to access the mass hierarchy, as long as the baseline is 
long enough. In this case, Ax 2 > 1 can be achieved, and one would require 6.8Mtonyr of 
data for a 3a signal with the TASD. Thus, given a detector ten times as big as NOvA, one 
could determine the neutrino mass hierarchy via z/ M disappearance after 23 years of data 
taking, but without pushing the technical limits to the extreme. 

One could also use a NOz/A-like superbeam experimental setup [41] and obtain similar 
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Figure 8: A% 2 for the test of the normal hierarchy as function of the simulated Am^ 
and baseline L of the SB-Upgrade, where the SB-Upgrade is operated with two years of 
neutrino running and six years of antineutrino running and not combined with any other 
experiment. The vertical lines correspond to the current best-fit value (right) and lower 
end of the allowed range (left, 90% CL [35]). The numbers give the A% 2 for the respective 
contours, where it is increased in steps of 0.5 (except first contour). 

quantitative results. Since the NOi/A superbeam has an energy spectrum peaking around 
~ 2 GeV, one can read off from Fig. |5] (right) that this experiment could provide the nec- 
essary information about the higher oscillation maxima, too. For example, the required 
luminosity for the NuMI beam with the Fermilab proton driver (see Refs. [11, 12] for the 
physics potential) would be around 19Mtonyr (for 3cr) using a comparable detector tech- 
nology to the neutrino factory in the same off-axis configuration as NOz/A. In this case, the 
optimal baseline would be around 5 500 km (similar to the neutrino factory), and neither the 
performance nor the baseline optimization would be significantly affected by a better energy 
resolution. Because of the conceptual similarity to the neutrino factory, we do not follow 
this approach anymore, and use "NB-Superbeam" to refer to a low energy (E < 1 GeV) 
narrow band beam. 

Though we have found that a lot of luminosity is required, the final increment required after 
a long fruitless search for nonzero sin 2 (26*13) effects may not be completely unrealistic. As 
depicted in Fig. [JJ there are many potential detector locations for such a long baseline and 
different potential neutrino factory laboratories. For example, detectors could be located in 
the Soudan mine for neutrinos coming from CERN or RAL, in India from JHF or CERN, 
in Finland from FNAL or BNL, in the Gran Sasso laboratory from BNL, or in Australia 
from JHF. 

We now discuss the robustness of these results. As we will see in the next section, the energy 
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resolution of the superbeam detector can drastically change its potential to determine the 
mass hierarchy, as well as the optimum baseline for doing the job. Hence, as already 
advertised, the detector performance is the "bottleneck." For the neutrino factory, on the 
other hand, it is the small number of events at low energies that limits the performance: 
One can see in Fig. El (right) that a factor of two better or worse energy resolution does not 
change qualitatively the potential to resolve the highest significantly contributing oscillation 
maximum. For AE > 250 MeV, however, we find a substantial degradation of the mass 
hierarchy sensitivity (larger than 50% change in Ax 2 ). In addition, the energy resolution 
hardly affects the baseline optimization of the neutrino factory. 

Another important issue is the precise knowledge of the oscillation parameters other than 
Am^. We find that perfect knowledge of these parameters would translate into a marginal 
20-30% reduction the luminosity requirements for all setups considered above. Finally, the 
impact of the hierarchy (sensitivity to normal or inverted hierarchy) turns out to be rather 
small and not to affect the optimization discussion. In summary, the energy resolution is the 
limiting factor for the NB-Superbeam, while the lack of statistics at low energies limits the 
ability of the neutrino factories to uncover the neutrino mass ordering in the limit #13 — > 0. 

Before proceeding, it is interesting to comment on the dependence of the sensitivity to 
the mass hierarchy on the simulated value of Am^. Fig. |H] depicts the A% 2 for the test 
of the normal hierarchy as function of the simulated Am^ and baseline L of the NB- 
Superbeam. In this case, we choose an option with two years of neutrino running and 
six years of antineutrino running not combined with any other experiment (one experiment 
only). As discussed above, the A% 2 at the current best-fit value of Am^ requires substantial 
luminosity upgrades. If, however, Am| x was at the lower end of the currently allowed 
region, a factor of four higher A% 2 can be obtained at a somewhat longer baseline. For 
even smaller values of Am| 1; even a 90% CL hint would be possible with this experiment 
only, where the optimal baseline increases with decreasing Arn^j . The reason for this better 
performance is that the solar and atmospheric mass-squared differences are getting closer for 
smaller values of jAm^l, and thus any effect which requires the presence of both oscillations 
becomes enhanced. The longer baseline becomes necessary because of the fixed peak of 
the energy spectrum where oscillations can be resolved with the chosen energy resolution 
the relationship Am^ x L ~ const, (weighted with 1/L 2 event rate drop) applies (cf., 
Eq. ©). On the other end of Fig. |Hl we find a substantial loss of sensitivity above Am^ ^> 
2.5 • 10~ 3 eV 2 , which means that the mass hierarchy measurement with a NB-Superbeam 
becomes very hard if Am| x turns out to be much larger than the current best-fit value. 
In this case only a neutrino factory approach could be possible, as it turns out that the 
neutrino factory performance and optimization is rather unaffected by the simulated value 
of Am^ in the currently allowed Amf^-range (90% CL). 

5 Pushing the limits 

From the last section, it is clear that the neutrino factory capability to determine the mass 
hierarchy in the limit #13 — > is hard to improve because of the limited statistics at low 
energies. One could, of course, always think about making use of an even lower energy 
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Two NB-Superbeams, v's only 



Two NB-Superbeams, v (2yr)+v (6yr) 



NB-Superbeam+NuFact 




Figure 9: A% 2 for the test of the normal hierarchy as function of two baselines for the 
different experiment combinations as given in the plot labels. In this figure, the energy 
resolution of the SB- Upgrade is pushed by a factor of two beyond the current understanding. 
For the oscillation parameters, the values in Eq. (j^J) are used, and Am^ = +2.2 10 -3 eV 2 . 

neutrino factory. We, however, regard this as qualitatively similar to a narrow band beam 
peaking at lower energies and therefore choose the NB-Superbeam as a representative for 
a beam with good low-energy statistics. In this section, we discuss the consequences of 
pushing the detector technology beyond "current" limitations. 

Fig. El depicts the combination of two NB-Superbeams (left - eight years neutrino running 
only), two NB-Superbeams (middle - two years v running and six years v running), and a 
NB-Superbeam with the neutrino factory, as defined in the previous section (right). In this 
figure, the maximum energy resolution of the NB-Superbeam(s) was reduced by a factor 
of two with respect to the Fermi motion-dominated effective energy resolution of about 
85 MeV. Comparing Fig. El to Fig. El we can read off that the baseline optimization for 
the longer NB-Superbeam baseline changes. In this case, the optimum is at L ~ 4 000 km, 
longer than the value obtained with worse energy resolution depicted in Fig. El The reason 
for the increase of the optimal baseline with the improved energy resolution, clearly observed 
here, was discussed in Sec. El 

Instead of using two NB-Superbeams at different baselines, a NB-Superbeam experiment 
with a single baseline which takes advantaged of combined neutrino and antineutrino run- 
ning can exceed the 90%CL threshold, as depicted in Fig. El (middle). Similarly, the com- 
bination with a neutrino factory instead of a second superbeam would yield even better 
results because of the high statistics of the neutrino factory (see Fig. El (right)). Note that 
Fig. El (right) is "orthogonal" to Fig. El (middle), which means that, here, the neutrino fac- 
tory plays a secondary role when it comes to determining the mass hierarchy (even though 
one NB-Superbeam alone is not capable of providing a 90%CL signal). By comparing the 
left and right panels of Fig. El we find that, say, a shorter baseline L ~ 3 000 km neutrino 
factory (plus the appropriate detector) could replace the role of the shorter NB-Superbeam 
baseline in order to yield similar results. Such a setup may already be available by the 
time the measurement discussed here is contemplated. In summary, we learn from Fig. El 
that an improved energy resolution would make the NB-Superbeam the key experiment for 
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Figure 10: A\ 2 for the test of the normal hierarchy as function of the energy resolution 
AE and baseline L of the SB-Upgrade, where the SB-Upgrade is combined with a neutrino 
factory at 6 000 km. The thick curve represents the optimum baseline for each individual 
energy resolution. The vertical line corresponds to the Water Cherenkov detector's effective 
energy resolution dominated by Fermi motion. For the oscillation parameters, the values in 
Eq. © are used, and Am^ = +2.2 10" 3 eV 2 . 

the mass hierarchy measurement for sin 2 (2^13) = 0, where the optimal long baseline hardly 
depends on the combination of experiments considered, but depends more strongly on the 
achievable energy resolution. 

How good an energy resolution do we really need in order to significantly boost the capa- 
bilities of a NB-Superbeam to determine the neutrino mass hierarchy in the limit #13 — > 0? 
Fig. El depicts A\ 2 as function of the energy resolution AE and the baseline L of the 
SB-Upgrade, where the SB-Upgrade is combined with a neutrino factory at 6 000 km. From 
this figure, one can read off a number of different results: First, the optimal baseline of the 
NB-Superbeam becomes longer as the energy resolution improves (thick curve). Second, 
for the energy resolution of the Water Cherenkov detector (vertical line), there is no sta- 
tistically significant signal and tremendous luminosity upgrades are required. However, for 
better energy resolution, the statistical significance is markedly improved. A factor of four 
better energy resolution would allow one to obtain a 3a determination of the mass hierar- 
chy for realistic luminosity levels. And third, the gradient in the direction of the energy 
resolution is rather significant, which means that any improvement of the energy resolution 
would reduce the required luminosity for a 3a signal. Note that though Fig. [TDlis computed 
for the combination of the NB-Superbeam with the neutrino factory, these results are qual- 
itatively unaffected by the particular combination of experiments (assuming it includes a 
NB-Superbeam) . 
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How would one go about improving the energy resolution? A detailed study of this issue 
is beyond the intentions of this paper, but we wish to add a few qualitative remarks. As 
already alluded to above, the source of energy resolution boils down, at the levels in which 
we are interested, to the Fermi motion of the nucleons inside the nuclei that make up the 
target. In a water Cherenkov detector, for example, sub-GeV neutrinos scatter, mostly 
quasi-elastically, off of nucleons contained in oxygen nuclei. The Fermi motion of these 
leads to the intrinsic AE ~ 85 MeV quoted above. There are, at least in principle, a 
few ways to improve on this. Different materials have different intrinsic uncertainties on 
the nucleon energies. One naively expects that materials with larger mass number (A) 
would have smaller intrinsic AE values. One the other hand, detectors that are capable of 
fully reconstructing the final state hadronic state in an event-by-event basis (in the case of 
quasi-elastic events, fully determining the energy and momentum of the final state daughter 
nucleon) may bypass this uncertainty. Furthermore, if one could afford a pure proton 
(or even electron!) target, Fermi motion would be absent (replaced, perhaps, by much 
less significant thermal motion) and one could hope to improve on the energy resolution. 
Finally, energy resolution issues could be effectively erased if mono-energetic muon-neutrino 
beams were available (such as neutrino beams from muon capture on nuclei, similar to [46], 
or neutrino beams from stopped pion decay). In this case, one could repeat the analysis 
done here by placing different experiments at several different baselines in order to measure 
Am^ as a function of L/E. 

While some of these propositions sound intriguing, we remind readers that in order to per- 
form the measurement we are proposing, extreme luminosities are required: not only do 
we need a more precise detector technology, but we also need one that is capable of accu- 
mulating data as quickly as the "bread-and-butter" designs which are currently available. 
Our hope here is not to resolve this issue in any way, but to attract the attention of the 
community to these matters. 

6 Summary and conclusions 

Determining the neutrino mass hierarchy is of the utmost importance, but may prove very 
challenging if 13 is too small. In the near future, we hope to learn, through next-generation 
long-baseline and reactor experiments, whether sin 2 2#i3 > 10~ 2 (see, for instance, Ref. [47]). 
If this is the case, these (or upgrades) should be able to determine the neutrino mass hierar- 
chy. Otherwise, we will need to wait for "next-next" generation probes, which are sensitive 
to sin 2 (26 , i 3 ) down to sin 2 20 13 ~ 10 -4 by, for instance, using the "magic baseline" at a neu- 
trino factory [18]. Again, if this turns out to be the case, these (or, say, a nearby supernova 
explosion) should be able to determine the neutrino mass hierarchy. If, on the other hand, 
#13 is smaller still, we are "stuck" with detailed studies of muon neutrino disappearance [30], 
and, perhaps, "non-oscillation" probes of neutrino masses (see, for example, [48-50]). 

In this study, we have re-visited the requirements for a neutrino mass hierarchy measure- 
ment using P m , in the case sin 2 (26*i 3 ) = 0. We need two very precise measurements of the 
atmospheric oscillation frequency at very different values of L/E, one at a very high oscilla- 
tion maximum where solar effects contribute significantly (alternatively, a combination with 
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antineutrino running could replace the combination with a precise Am^ measurement from 
a "small" L/E value). Theoretical estimates (analysis at the probability level) point to the 
lowest attainable neutrino energies and very long baselines, L > 3000 km. In practice, on 
the other hand, the detector energy resolution limits our ability to resolve fast oscillations 
(and hence "measure" Am^) and the optimal baseline is determined to be such that the 
largest number of events is obtained at relevant oscillation maxima. The conditions above 
require large event rates capable of overcoming the large 1/L 2 suppression of the beam 
flux with which we are required to deal, and excellent energy resolution in order to resolve 
the relevant low-energy oscillation maxima. We were especially interested in the whether 
"current" experimental setups meet these requirements, at least in principle. Not surpris- 
ingly, we find that, in order to significantly determine the character of the neutrino mass 
hierarchy, one is forced to push existing detector capabilities (and running times) to their 
technological extremes. 

Using conventionally established detector technology, we have demonstrated that a neu- 
trino factory optimized for this purpose, operated at L ~ 6 000 km and with substantially 
increased statistics could do the mass hierarchy measurement for sin 2 (2#i 3 ) = 0. From the 
optimization, we find it best not to require charge identification in order to enhance the 
low energy detection efficiencies, and we use a low muon energy = 15 GeV to relatively 
enhance the flux in the low-energy region. For example, a TASD detector ten times as big 
as NOvA would allow, after 23 years of running time, a 3a measurement (for 10 21 useful 
muon decays/year). Though these numbers are quite large, note that the mass hierarchy 
measurement for sin 2 (2#i3) = becomes most relevant if no sin 2 (2#i3)-signal is found. As- 
suming that substantial resources will go into this (assumed to be fruitless) search for 9iz, 
it is reasonable to suppose that a significant amount of data and equipment (beams and 
detectors) will already be available at the time we decide to measure the mass hierarchy as 
described here. In addition, the 6 000 km baseline may also be required for a high-confidence 
level MSW effect verification for sin 2 (2#i3) = via the solar appearance term [51]. 8 The 
neutrino factory measurement is limited by the flux-related drop of events in the low energy 
region, which is quite common for similar broad band beams. This means that in the region 
where the energy resolution limits the performance, the event rates are already too low to 
have statistical significance, and a better energy resolution would not help significantly. 

A very promising alternative is a narrow band off-axis superbeam similar to T2HK. Using a 
conventional water Cherenkov detector, the effective energy resolution coming from Fermi 
motion of the nucleons limits the ability to resolve higher oscillation maxima, and yields a 
performance much worse than that of the discussed neutrino factory. However, with modest 
improvements to the energy resolution, such a narrow band beam measurement combined 
with a measurement at a smaller L/E (or antineutrino running or a neutrino factory) proves 
to be significantly more sensitive to the neutrino mass hierarchy. For example, if the energy 
resolution were improved by only a factor of two (from about 85 MeV to about 43 MeV) 
one would be able to determine the neutrino mass hierarchy at the 90% CL by combining 
eight years of data from a T2HK-like experiment at L = 4 000 km (1 Mton Water Cherenkov 

8 Notc that for this appearance channel measurement a detector with charge identification is required. 
Therefore, one may want to a priori choose a detector technology which can be magnetized at this baseline 
and analyze the data set in two different samples (with and without charge identification). 



18 



detector) with data from a neutrino factory at L — 6 000 km (50 kt TASD). We note that 
the energy resolution limits the number of the observable oscillation maxima and that it 
determines the optimum baseline for performing the measurements described here. For a 
fixed low-energy narrow band beam (peak energy below 1 GeV), we find that the optimal 
baseline increases (up to around 6 000 km) as the energy resolution improves. 

In summary, we have identified two important discriminators between the use of low-energy 
narrow band beam technology and a relatively low energy neutrino factory: If Am| x turned 
out to be substantially below the current best-fit value, or a detection technology with a 
better minimum energy resolution than that of Water Cherenkov detectors was identified, 
narrow band superbeams should be the appropriate choice for a #13 = mass hierarchy 
measurement. If, on the other hand, Am^ was larger than the current best-fit value 
or no substantial improvement of the energy resolution was obtained, a ("lower-energy," 
= 15 GeV) neutrino factory at about 6 000 km with a substantial luminosity increase, 
or another beam with similar characteristics, could be the most efficient choice for the mass 
hierarchy measurement. 

We conclude by reemphasizing that energy resolutions not far from the currently achievable 
would significantly enhance the physics potential of this measurement. For example, differ- 
ent detector technologies and/or detection materials (such as liquid argon) could be used to 
improve the energy resolution and our ability to determine the neutrino mass hierarchy - 
provided that the ability to record events is comparable to the QE event rate in a 1 Mton 
Water Cherenkov detector. Since improving on the low-energy muon energy resolution has, 
so far, not been necessary for the more conventional measurements, we wish to raise the 
level awareness of the community to this issue. 
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